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In models like axion monodromy, temporal features during inflation which are not associated with 
its ending can produce scalar, and to a lesser extent, tensor power spectra where deviations from 
scale-free power law spectra can be as large as the deviations from scale invariance itself. Here the 
standard slow-roll approach breaks down since its parameters evolve on an efolding scale AN much 
smaller than the efolds to the end of inflation. Using the generalized slow-roll approach, we show that 
the expansion of observables in a hierarchy of potential or Hubble evolution parameters comes from 
a Taylor expansion of the features around an evaluation point that can be optimized. Optimization 
of the leading-order expression provides a sufficiently accurate approximation for current data as 
long as the power spectrum can be described over the well-observed few efolds by the local tilt and 
running. Standard second-order approaches, often used in the literature, ironically are worse than 
leading-order approaches due to inconsistent evaluation of observables. We develop a new optimized 
next-order approach which predicts observables to 10“ 3 even for AN ~ 1 where all parameters in 
the infinite hierarchy are of comparable magnitude. For models with AN -C 1, the generalized 
slow-roll approach provides integral expressions that are accurate to second order in the deviation 
from scale invariance. Their evaluation in the monodromy model provides highly accurate explicit 
relations between the running oscillation amplitude, frequency and phase in the curvature spectrum 
and parameters of the potential. 


I. INTRODUCTION 

In a general inflationary model, temporal scales that 
are not directly associated with the end of inflation leave 
their imprint in cosmological observables. For canonical 
single field models, these arise from features in the po¬ 
tential which leave the curvature and gravitational wave 
power spectra nearly scale invariant but no longer scale- 
free power laws. Current observations constrain devi¬ 
ations from scale-free power law spectra only as well as 
deviations from scale invariance itself [1] . Therefore these 
constraints only impact models with large features. 

In the standard slow-roll approximation, the leading- 
order effect would be a running of the tilt of which the 
amplitude is quadratic in the small deviation from scale 
invariance. It is common to test for a finite value of such a 
parameter in data sets like the Planck CMB power spec¬ 
tra [1]. Yet the running of the tilt is only constrained 
at the ~ 1CT 2 level, the same level as the tilt from scale 
invariance itself. Interpreting these constraints requires 
going beyond the standard approximation where slow- 
roll parameters are taken to be both nearly constant 
and strongly hierarchical (e.g. [2-5]). Using the stan¬ 
dard second-order approach on such data, as is common 
in the literature (e.g. [1, 6-9]), can provide misleading 
results for models where the running of the tilt is not 
exactly constant or the potential purely cubic. 

On the other hand, if the observed power spectra are 
well characterized by the local tilt and running of the 
tilt, which are the first two terms in a Taylor expansion 
of a continuous feature, then a slow-roll hierarchy still 
exists and can be used to predict observables accurately 
(cf. [10]). It is simply that the efolding scale of the feature 
AN, while still greater than unity, is much less than the 


number of efolds to the end of inflation. 

In this paper, we use the generalized slow-roll (GSR) 
approach [10-16] to extend the validity of the standard 
slow-roll approximation. Here the deviations from scale 
invariance are only assumed to be small in amplitude 
not in temporal frequency. With GSR, we explore the 
relationship between features in the potential and the 
hierarchy of slow-roll parameters. Based on a temporal 
Taylor expansion, we develop optimized approaches for 
the evaluation of curvature and gravitational wave power 
spectra from inflation. We test these approaches in the 
axion monodromy model [17] where low frequency cases 
produce a slowly varying running of the tilt [18-20] and 
high frequency cases imprint oscillations in the power 
spectra [21, 22]. 

The paper is organized as follows. In §11, we review 
the GSR approach and iterate it to second order in the 
amplitude of deviations from scale invariance to describe 
power spectra in terms of integrals over the temporal his¬ 
tory of inflation. We consider first order terms in §111 and 
show how the Taylor expansion of the temporal history 
is related to observables and to the hierarchy of slow-roll 
parameters of the potential and Hubble evolution. Op¬ 
timized evaluation of the leading-order Taylor term [11] 
provides predictions and inflation potential reconstruc¬ 
tions that are accurate to next-to-leading order. These 
suffice for models with running of the tilt that is nearly 
constant across the observable scales. In §IV, we gen¬ 
eralize the optimized approach to second order and an 
optimized next-to-leading-order approximation and con¬ 
trast these approaches with the standard second-order 
approximation which inconsistently treats these evalu¬ 
ations. We illustrate these optimized evaluation tech¬ 
niques for monodromy in both the low and high frequency 
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regime in §V. We discuss these results in §VI and pro¬ 
vide relationships between the various parameterizations 
in the Appendix. 


II. GENERALIZED SLOW ROLL 


The GSR approach is ideal for studying and extend¬ 
ing the validity of the standard slow-roll approximation 
as it assumes only that the amplitude of deviations from 
a de Sitter expansion are small with no assumptions on 
their temporal frequency or equivalently on the smooth¬ 
ness of the potential [10, 11, 13, 14]. In this approach, 
we iteratively solve the exact Mukhanov-Sasaki equation 
for the evolution of the inflaton mode function or field 
fluctuation in the spatially flat gauge 6<f> = y/av2k 


d 2 y 

dx 2 



y_ 

V / ) * 2 ’ 


( 1 ) 


to obtain the comoving curvature power spectrum 


A 


2 

n 


lim 

x —>-0 
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Here ' = d/d. In r/, x = krj and rj is the (positive, decreas¬ 
ing) conformal time to the end of inflation and the de¬ 
viations from de Sitter expansion are given by variations 
in 

f 2 = 8 ? T_tH_ , ( 3 ) 

where the Hubble slow-roll parameter is defined by 


dhiH _ 1 ( d(j)\ 2 

dN ~ 2 \dN J ’ 


( 4 ) 


with N as the (negative, increasing) number of efolds 
relative to the end of inflation. 

The fundamental assumption in this approach is that 
the mode function y remains close to its de Sitter form 


y 0 (x) = fl + £) e ix , (5) 

and we will loosely refer to this property as requiring the 
amplitude of deviations from de Sitter or scale invari¬ 
ance encoded in / to be small (see [23-26] for exceptional 
cases). If so, we can take the formal Green function so¬ 
lution to Eq. (1), 

y(x) = yo(x)- f ^-— -^-y(u)lm[yZ(u)y 0 (x)], (6) 

replace y —>• yo on the right-hand side and iteratively 
improve the solution. 

To first order in the de Sitter deviations, this procedure 
results in [27] 

r°° dr 

lnA^j(fc) ~ / 0 = G(lnai*) + / — W(x)G'(lnx), (7) 

J cr .... *■£ 


where x* < 1, 


G = -21n/+?(ln/) / , 


( 8 ) 


and 


W{x) 


3sin(2x) 

2a; 3 


3cos(2a;) 3sin(2x) 
x 2 2x 


( 9 ) 


While we assume a canonical scalar field throughout, all 
of the GSR-based results here and below can be read¬ 
ily generalized to P(X, <j>) theories and the effective field 
theory that parameterizes them [16, 28]. 

The function W determines exactly how inflaton fluc¬ 
tuations freeze into curvature fluctuations as they pass 
the horizon. We can make the connection more explicit 
by integrating Eq. (7) by parts to obtain [27] 


r°° dx 

J 0 » - / —W'{x)G(\nx). (10) 

Jo X 


Note that G is a function of time In ?y alone and denot¬ 
ing it as a function of In x is simply a convenient choice 
of its zero point for a given k. Scale dependence of the 
power spectrum arises only if G' 0 and hence its value 
quantifies the deviations from de Sitter results. These 
typically must be 


G' = c(i) (11) 

in order for features not to prematurely end inflation. 
The overbar here represents an average over several efolds 
and denotes the fact that only the integral of G' needs 
to be small or equivalently that transient effects do not 
necessarily end inflation. This average as we shall see is 
closely related to the average tilt of the curvature spec¬ 
trum. In fact Eq. (10) is valid even if the inflaton poten¬ 
tial contains discontinuities or delta function sources in 
G' where the local tilt and average tilt differ substantially. 
Given that linx^o W' = 0, this form also conserves cur¬ 
vature fluctuations outside the horizon and gives a man¬ 
ifestly positive definite power spectrum, and it remains a 
controlled approximation for up to order unity deviations 
[14] unlike related variants [11, 13]. 

In particular, linearity in G' should not be conflated 
with the constancy of G'. If we define AN as the typical 
efolding scale of its variation, then 

= 0 (av) g ' = ° (iviv) < 12 » 

rather than 1/N 2 as the standard slow-roll approxima¬ 
tion assumes by requiring that the only temporal feature 
during inflation be associated with its end. We shall show 
in §111 that if 1 < AN <C |JV|, a slow-roll hierarchy of 
successively smaller derivative parameters still exists, but 
order counting needs to be generalized using a consistent 
Taylor expansion of G between observables. If AN <C 1, 
the hierarchy is inverted so that computing observables 
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requires resumming an infinite series or equivalently the 
direct integration of Eq. (10). In this regime observ¬ 
ables depend on the continuous function G (see §VB), 
and model-independent approaches seek to reconstruct 
that function rather than measure a series of parameters 
[29-31]. 

These G(l/NAN) effects should be distinguished from 
the true 1/N 2 ones. In the GSR approximation, these 
come from terms that are quadratic in the deviation from 
de Sitter solutions G'. Iterating the Green function ap¬ 
proach, we obtain [13, 14] 


A 


2 

n 


O lo 


1 + -T-tf 




(13) 


where the first-order term Iq is defined in Eq. (7) and the 
second-order corrections are 


1, = 


A. 


—G' (In x) X (x), 
x 


1-2 = —4 


dx 


1 / 


duf 


—\ x +i x )jJ -,J, ( 14 ) 


with 


X(x) = —z-(sina; — a: cos a;) 2 . 


(15) 


In the following sections we keep order counting in 1/N 
and 1/AIV distinct. 

For tensor fluctuations, the same GSR approach holds 
with the replacement of / with 


fh = 


(16) 


in the construction of the source Gh of the tensor fluctu¬ 
ations [15, 16]. With these substitutions, Eq. (13) then 
provides the tensor power spectrum in each polarization 
state In A^_ x (fc). 


III. OPTIMIZED LEADING- AND 
FIRST-ORDER APPROXIMATION 


A. Smoothness Hierarchy 


The hierarchical structure of the slow-roll parameters 
stems from a smoothness assumption for the deviations 
from a de Sitter expansion. For any smooth source of 
deviations G(lnx) in the GSR formalism, we can Taylor 
expand its form around an epoch near horizon crossing 
lmcy and integrate its effect term by term. The standard 
slow-roll approximation proceeds by assuming lnx/ = 0 
but we shall see there are advantages to tuning this evalu¬ 
ation point to make it correspond better to the freeze-out 
epoch. Note that a shift in this point corresponds to a 
shift in efolds of SN ss —5 lnx/. 

For the first-order lnA^ ss Iq term of Eq. (10), this 
leads to 

OO 

In A^ ps G(lnx/) + ^ g p (lna;/)G^ p ^(lnx/), (17) 

p=i 


where G ^ denotes the pth derivative of G with respect 
to In a; and 


QpQnxf) = -— [ ^W'(x) fin Jtl] . (18) 

p\ Jo x V x f J 


These coefficients can be calculated using the generating 
function [11] 


F(z,x f ) 



(2 xj) z cos 



3T(2 + z) 

(1 — * 0 ( 3 - 2 )’ 


(19) 


so that 


q p (lnx f ) 


i lim a-rc,*,) 

p\ z-x) dz p 


( 20 ) 


To clarify the dependence of these coefficients on lmry, 
it is also useful to express q p (\nxf) in terms of the first 
term 


qi(lnxf) = lnaq — lna;/, (21) 


Utilizing and extending the techniques of Ref. [11], we 
elucidate the conditions under which a slow-roll hierarchy 
of parameters for scalar and tensor power spectra observ¬ 
ables exists and the relative size of terms in the series in 
§IIIA. In §IIIB, we show how the next-to-leading-order 
terms in the hierarchical expansion are generally large 
whenever running of the tilt is comparable to the devia¬ 
tions in the tilt itself but can be absorbed into an opti¬ 
mization in the time of fluctuation freeze-out consistently 
between observables. The resulting optimized leading- 
order description suffices for scalar and tensor spectra 
that can be described by a nearly constant running of 
the tilt even when it is of order the deviations from scale 
invariance represented by the tilt itself (cf. [10]). We re¬ 
late this hierarchy to the Hubble slow roll parameters in 
§IIIC and to the potential slow-roll parameters in §IIID. 


where 

lnaq = ^ - In2 — r y E . (22) 

Here ")e is the Euler-Mascheroni constant. Since 


dq p 

dqi 


dq p 

dlnxf 


Qp~ 1 5 


(23) 


q p is a pth degree polynomial of q-[. Each degree intro¬ 
duces an extra constant of integration that is indepen¬ 
dent of In Xf. It is convenient to define these in terms 
of g p (lna:i) since gi(lna;i) = 0. The higher coefficients 
become 


9p(lna:/) 


Qi( l ^Xf) 

p\ 


+ g fc qh 1£l l 9 „ (lni;) 

n —0 


(24) 
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The next two terms are thus given explicitly by 


42 (In x f ) 
43(lili/) 


4i(hi x f ) 
2 

ViO-Wf) 

6 


4 — 3n 2 
72 

4 — 37r 2 
72~ 


4i (I 11 a;/) 


55 

81 


(25) 

C(3) 
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We optimize the slow-roll approximation below by choos¬ 
ing In Xf to set certain coefficients to zero. Equation (23) 
implies that setting q p (\nxf) = 0 makes q p +i(lnxf) take 
on extremal values as a function of In Xf. Certain opti¬ 
mizations also have the benefit that the extremum is a 
minimum of \q p +\ | and hence suppress the next correc¬ 
tion. 

Since 


dG^ (In x f ) 
dink 


-G^ +1 \lnx f ), 


the tilt and the running of the tilt are given by 


(26) 


d In A)/ 
din 


-G'(lnx f ) - ^ qpG^+V (In x f ), 

P= 1 


dn s 

a =- 

din A: 


G"(lnx f ) + J2q p G {p+2) {lnx f ), (27) 
p= 1 


with the obvious continuation to the running of each suc¬ 
cessive quantity. Since the q p coefficients are fixed given 
a choice of In Xf for all lnfc, we omit the lncc/ argument 
of q p for clarity and where no confusion will arise we 
also do so below for compactness. A nearly constant tilt 
in In A; requires \G"\ <C |G'|, a nearly constant running 
G'" -C \G"\, and each gains its next-to-leading-order 
correction from the q\ term. 

For tensor fluctuations, we similarly have [15, 16] 

r 00 dr 

In A 2 +<x (k)*~ —W'(x)G h Qnx) 

Jo x 

OO 

~ G h (lnx f ) + ^2 qpG^ (In x f ), (28) 
p =1 


which likewise determines the tensor tilt and running of 
the tilt 


n t = 


dlnA^_ x 

d In k 

dn t 


ot t = 


d In k 


-G' h (lnxf) - ^2q p G { ^ +1) (Inx f ), 

P= 1 

OO 

G^lnxf) + '£q p G ( * +2 \lnx f ).(29) 

p=i 


For a sufficiently smooth G(lnx) and G^(lna;), the se¬ 
ries expansion of the power spectra will rapidly converge. 
Since 


lim 

P->oo q p _i 


1 

2’ 


(30) 


the criteria for convergence for the scalar power spectrum 
is 



LO/SO 

OLO 

ONO 

GSR 

p 

0 

1 

2 

(X) 

AN 

> 50 

> few 

> 1 

all 

lna;p 

0 

1.06 

0.22 

- 

4i 

1.06 

0 

0.84 

- 

42 

0.21 

-0.36 

0 

- 

43 

0.10 

0.28 

0.078 

- 


TABLE I. Slow-roll approximations are characterized by their 
order p in a temporal Taylor expansion which determines their 
applicability in describing features spanning AN efolds. Stan¬ 
dard second-order (SO) approaches only improve on leading 
order (LO) by keeping p > 0 in some but not all observ¬ 
ables. Optimized evaluation (OLO/ONO) achieves consistent 
pth order accuracy with only p — 1 additional parameters by 
evaluating them at special values of In a:/ = lnx p . The coef¬ 
ficients q n for n > p control the error from truncation. The 
GSR approximation forgoes the Taylor expansion for an inte¬ 
gral approach. 


and is similar for the tensor spectrum and tensor source 
G;,. The series is dominated by the first term for 
sufficiently small values for this ratio, with each suc¬ 
cessive term suppressed by the smoothness scale AN, 
|Gbd/G( p-1 )| = 0(AN~ 1 ). Therefore, the Taylor ex¬ 
panded form above is applicable for features with AN > 
1. For high frequency features where AN <1/2 and the 
convergence criteria are violated, the GSR integral for¬ 
mula must be developed on a case by case basis (see §V B 
for the monodromy example). 


B. Optimized Leading vs. Next Order 

We call an approximation scheme that just retains the 
GSR Jo term and no corrections from the Taylor series, 
first order in deviations from scale invariance and leading 
order (LO) in the hierarchy. Namely, 

hi A^ « G(lna;/), 
n s - 1 « -G'(lna:/), 

a ~ G"(Inxf), (LO) (32) 

and similarly for the tensor observables. One that adds 
p = 1 terms we call next(-to-leading)-order (NO) in the 
hierarchy. These distinctions still allow us to choose the 
evaluation epoch In x / which can be exploited to make a 
specific LO approximation as accurate as the generic NO 
approximation. 

In the standard slow-roll approximation, one takes 
lna;/ = lna;o = 0 and to NO 

lnA^ m G(0) + 5i( 0)G , (0), 
n s —1«—G'(0)- 9 i (0)G"(0), 

a«G"(0) + (?1 (0)G , "(0). (33) 


lim 

p—too 


G(p) 


< 2 , 


(31) 


Here we have restored the In x / argument of q n for clarity. 
Since gi(0) ~ 1.06 (see Table I), if a is comparable to 
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n s — 1 then the NO corrections are comparable to LO 
and are therefore required for accuracy. 

On the other hand, since all observables follow this 
same Taylor series form, this NO approximation is equiv¬ 
alent to shifting the evaluation epoch by <71 (0). This 
brings the evaluation point to 

In 27 = gi(0) « 1.06, (OLO) (34) 

which using Eq. (21) is equivalent to setting 

q^lnxf) = gi(lnari) = 0. (35) 

We call this the optimized leading order (OLO) approx¬ 
imation. By adopting this optimization, we gain all of 
the benefits of a next order approximation without any 
of the complexity. From this point forward, unless we 
specify otherwise, the LO approximation will refer to the 
standard LO approximation of In Xf = 0. 

Beyond the NO approximation of Eq. (33), this shift 
does not resum the terms in the Taylor series since for 
p > 1 

<fr(0) ^ ® (36) 

or equivalently for the OLO evaluation point < 7 p (lnxi) 7 ^ 
0. Furthermore |<72 (In ) | is a local maximum of 

|g 2 (Ina-’/)I so that zeroing the NO correction comes at 

the expense of maximizing the next-to-NO correction. 
For a low frequency AN 1, i.e. for a strong hierar¬ 
chy \G^/G^ P ~ 1 ' > | <C 1, this is a small price to pay. For 
AN ~ 1, high accuracy requires going beyond the OLO 
optimization. In §IV we will optimize the NO approxi¬ 
mation itself by zeroing the next-to-NO correction from 
<72 ■ We shall see that this has the added benefit that one 
can choose the solution that is a local minimum of | < 731 . 

It is important to note that this series of optimizations 
is not equivalent to truncating at a fixed order in 1/AAb 
Since 


= (37) 

dropping 0{\/N AN) would eliminate G" in all observ¬ 
ables including a and dropping 0(1/N AN 2 * ) would re¬ 
tain the q± correction for n s — 1 but not a. By truncating 
at the same order in the Taylor expansion of each observ¬ 
able, we have ensured that they are all consistently eval¬ 
uated at the same effective epoch when the corrections 
are resummed. We shall see that this feature is crucial 
for ensuring consistency between observables. Table I 
summarizes the various approximations and their appli¬ 
cability. 

Since derivatives with respect to In 77 and In k are in¬ 
terchangeable, the accuracy of this truncation at NO or 
equivalently OLO is directly related to the accuracy with 
which the power spectrum can be described by the local 


tilt and running around a given pivot scale ko 
In A \{k) « In A s + [ n s (k 0 ) - 1] In 

+ ^V(A). ” ,38) 

If the observed power spectrum is a good fit to this form 
over at least the AN ~ few that are observationally well 
constrained, then OLO will also provide a good approx¬ 
imation. Omitted corrections are suppressed relative to 
OLO by 1/AiV 2 for each observable .* 1 We shall quantify 
this consideration in §V A. 


C. Hubble Slow Roll Parameters 


Although the G^(lnx/) terms can themselves be 
thought of as the hierarchy of slow-roll parameters that 
are the most directly related to observables, it is useful 
to reexpress them in terms of the more familiar Hubble 
slow-roll parameters tjj and 


1 (d 
~ HP(j) \dt 



(39) 


or equivalently derivatives of In H with respect to efolds 
through the hierarchy relations 


din en 
dN 
d In S p 

dN 


— 2 (en + Si), 


°p+ 1 


+ pe H - 6 i- 


(40) 


Using these relations we can express G < - p> in terms of the 
Hubble slow-roll parameters eu, S p as detailed in the Ap¬ 
pendix. In particular, the hierarchy of G^ derivatives 
is equivalent to the hierarchy of S p parameters. Equa¬ 
tion (40) says that a hierarchical scaling with |<5 p+ i/<5 p | < 
1 requires the fractional change in S p per efold to be small 
in addition to e#, |<5i | <C 1. In particular, the Taylor ex¬ 
pansion in G^ becomes a summation over the hierarchy 
S p (see the Appendix) 


G (P) « (-l) p (25p + ?<w) + O , (41) 

for p > 2. For order counting purposes 

«- 4 ' = °(v)’ 4 - = °Gv5V^)- < 42 > 


Although our expressions for observables in terms of 
G are explicitly 0(G') = G(N~ 1 ) due to the linearity 


1 The apparent contradiction with the power law example in 

Ref. [10] is due to numerical problems in their calculation [E. 

Stewart (private communication)]. Their feature model calcula¬ 

tion is also in error. 












6 


of Eq. (17), we can choose whether or not to keep non¬ 
linear terms in the conversion of Hubble slow-roll pa¬ 
rameters to G. Different choices will differ by 0{N ~ 2 ) 
and hence include a subset of second-order corrections 
arising from the evaluation of the first-order GSR term. 
Our convention will be to convert the OLO expression to 
0(N~ 1 ) in In and n s — 1 but retain 0(N~ 2 ) terms in 
a. This is to maintain backward compatibility with the 
standard slow-roll approximation where leading order for 
a is assumed to be 0{N~ 2 ). We also drop terms that are 
0{N~ 2 /AN). Thus, our leading-order expression spans 
the range of possibilities from AN -C \N\ to AN ~ \N\, 
i.e. large to small running of the tilt. In this sense, our 
OLO approach is more general than the optimization in¬ 
troduced in Ref. [11] where AN <C |1V| is assumed. 

With this convention the OLO scalar observables are 




10 2 
“ ~J €H ~ 3 51 


(43) 


n s - 1 « -4e ff - 2Si - -S 2 


(OLO) 


a ~ —2<$2 — — <$3 — 8 ejy — 10e#<5i + 2S 2 

O 


Here and below the notation | applies to the whole ex¬ 
pression. These differ from the standard slow-roll approx¬ 
imation in that <5 p+ i terms are not dropped just because 
a lower-order S p appears. 

Likewise the tensor observables become 


lnA^ jX 


n t 

a t 


, H 2 8 

n 2k 2 ~ ?> eH 


-2e H 


— 4 e 2 H - 4 :£ H Si 


(OLO). (44) 


These expressions for n t and at are the same as the stan¬ 
dard slow-roll expressions due to the lack of a running of 
n t at linear order in e# and 5 n . 

Since OLO self-consistently accounts for both running 
of the tilt and running of the running of the tilt by ab¬ 
sorbing the next correction into the evaluation point, it 
is both a simpler and better approximation than keep¬ 
ing all 0{1/N 2 ) and 0(1/N AN) terms in the standard 
slow-roll approximation as often used in the literature 
(e.g. [1, 6 , 7]; see §V for tests). Note also that the 
standard approach is often phrased in terms of the Hub¬ 
ble flow parameters where 5 P and its hierarchy equation, 
Eq. (40), is replaced by 


dlne p 

^7\T =ep+1 ’ 


(45) 


so that 


While the two series are algebraically equivalent when 
all terms are kept, Hubble flow parameters are inconve¬ 
nient for cases that lack a strong hierarchy [see Eq. (A17) 
and Fig. 14]. Equation (45) says that for e p+ i to remain 
small, the fractional change in e p over an efold must be 
small. Hence order counting in powers of e p automati¬ 
cally conflates AN and N as Eq. (46) illustrates. For 
example in the cases considered in §V where the varia¬ 
tions occur on AN <C 1, these parameters would have 
poles, whereas the S p do not. 


D. Potential Slow Roll Parameters 


We can likewise relate the G ^ hierarchy to the equiv¬ 
alent in derivatives of the potential. Using the exact 
background equations of motion for the inflaton and the 
expansion 

(3 + 6 1 )H 2 ^ = -VW, 

(3 - e H )H 2 = U, (47) 


we can relate the derivatives of the potential to the Hub¬ 
ble slow-roll parameters 


U = 

Vi = 



1 /( 2 ) 


V 


_ 9 (1 + V3 ) 2 

“ H (l~e H /3) 21 
€H — Si — S 2 /3 
1 — £h/ 3 


(48) 


Note that to lowest order Eq. (48) gives U « 2ch and 
Vi » - 6 2 /3. As shown in the Appendix, these 

imply G' ~ 3U — 2Vi. Likewise the series of higher-order 
derivatives with respect to efolds is replaced by higher- 
order derivatives of the potential with respect to the field 


V p = 



P_1 y(p+ 1 ) 
V 


(49) 


where = d p V/d(j) p . In the literature the first few 
parameters are also known as ey = U/2 , rjy = Vi, and 
= V 2 - Terms that are quadratic in U,V P are higher 
order in the deviation from a de Sitter expansion. How¬ 
ever, the relative size of the V p parameters depends on 
the smoothness of the potential V((j>) just as that of S p 
depends on the smoothness of H(N). 

We explicitly carry out these conversions in the Ap¬ 
pendix and outline the results relevant to the G{l/N) 
expansion here. Note that to leading order Eq. (47) pro¬ 
vides the condition for a friction dominated roll or at¬ 
tractor solution, 


d(t> _ 1 U« _ yw 
dN ~ ~3lP ~ V~ 


(50) 


El = € H , 

£2 = 2 (en + <5i), 

£263 = 4 ejj + 6e#(5i — 2<5 2 -I- 2S 2 - 


which also converts the two senses of derivatives 

d _ yW d 

(46) din rj V d<j>' 


( 51 ) 
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or 


dU _ dV p 
dlnrj ’ d In ?] 


V p+ i, 


and thus 


G (p) » —2V P + 0(N~ 2 ), 


(52) 


(53) 


for p > 2. The Taylor series in G^ has a one-to-one 
relationship to the same series in V p (see the Appendix 
for details). The convergence criteria for the Taylor ex¬ 
pansion in Eq. (31) likewise becomes a direct condition 
on the potential slow-roll parameters 


lim 

p—> OO 


G(p) 

GKp- 1 ) 


lim 

p—too 


V v 


Vp-! 


< 2 . 


(54) 


We can interpret this condition by approximating 


V p _ d(f> V^+V _ dliiV^ 

~ ~dN V(p) dN ' ( ^ 

Thus the convergence condition is that the fractional 
variation in V ^ experienced by the field across an efold 
is small. 

Although the OLO approximation again follows di¬ 
rectly from evaluating G^ in terms of these parameters, 
it is useful to restore the general evaluation epoch In a;/ 
to clarify its relationship to the field position on the po¬ 
tential. Equations (17) and (27) become 


A 2 ~ V 
11 12t x 2 U 


l+hq 1 ~-)u-2Y / QpV P 


P= 1 


n. 


-l«-3W + 2V 1 + 2^g p V p+1 , 


P= 1 


a : 


-2V 2 - 2 J2 QpVp+2 - 6U 2 +SUV!, (56) 

p=i 


where the potential terms are evaluated at the field po¬ 
sition 0(ln Xf) which we will write as In Xf as shorthand 
notation. In this language, optimization of Ins/ is equiv¬ 
alent to optimization of the field position at which to 
evaluate the potential parameters. From the field po¬ 
sition lira;/ = lnso = 0 , a shift in evaluation point to 
In si = <?i( 0 ) causes a shift of 


8<j> 


y(V yW 
-5N^— qi ( 0) 


V 


V 


(57) 


or a change in potential parameters 
5V k V^dfj) KUV qi (0), 

yW 

5U « 2 —(Vi - U)8(j> « 2W(Vi - U) q r(0), 

8V P » ^yV p+ i^ » V p+l9l (0) (58) 


that cancel the NO corrections in the hierarchy. In fact 
starting from an arbitrary lira;/ we can define 


V p (lns/) = V p (lns/) + g 1 (lns/)V p+ i(lns/) 

~ Vp(lnsr). (59) 


For the OLO choice, V p (lnsi) = V p (lnsi), and for al¬ 
ternate Ins/, V p (lns/) differs from Vp(lnxr) only by 
next-to-NO corrections. V n will be useful for compar¬ 
ing higher-order approximations in §IV. In fact for an 
arbitrary order we could define V p = V p + qn V n + P to 
absorb the appropriate number of correction terms into 
an effective potential parameter. 

Thus the OLO approximation simply amounts to shift¬ 
ing the evaluation point in field space by the same 
amount for all observables. Maintaining this consistency 
condition is important in testing models. Aside from this 
evaluation point, the OLO approximation (In Xf = lnsi) 
takes the same form as the standard LO approximation 

(In x f = 0 ), 


A 


2 

n 


n s - 1 
a 


1^7 

-3Z7 + 2 V,, 
-2V 2 -6W 2 + 8 WVi, 


(O/LO). (60) 


Though identical in form, the standard approximation 
assumes V 2 is 0(N~ 2 ) making a suppressed with respect 
to n s — 1 and absent in its leading-order prediction, not 
because of resummation to a new evaluation point but 
because of its truncation at 0(N~ 1 ). Again, the advan¬ 
tage of this OLO is that not only is it as simple as the 
standard LO slow-roll prescription, but it consistently in¬ 
corporates the NO corrections that run the tilt and run 
the running of the tilt. Corrections to it are suppressed 
by 0(l/AiV 2 ). 

We shall see that the same is not true for approaches 
that correct for evolution at a fixed order in powers of 
AN ~ N as is often used in the literature (e.g. [3, 8 , 9]). 
For instance if a is of order n s — 1, the standard second- 
order approach is really only appropriate for strictly con¬ 
stant a or equivalently a purely cubic potential where the 
cubic term is large but no higher terms are important. 
Importantly, in OLO we absorb the effect of the quar- 
tic term V 3 in the optimized evaluation of a to correct 
for the running of its value between horizon crossing and 
freeze-out. In the second-order approach, tilt is run to 
the freeze-out point, but running is not. One would then 
infer an incorrect relationship between the derivatives of 
the potential and hence potentially incorrectly falsify a 
true model that is not purely cubic from observations. 
Thus, a second-order approach is both more complicated 
and less general than the OLO approach. We illustrate 
this in §V A for the monodromy example. 

(p) 

For tensor fluctuations, conversion of G£ to the po- 





















A. General Expression 


tential parameters gives 


A 


2 

+,X 


(k) 


r 

n t 

at 


V 

6n 2 


X ~l U 


8 U, 

-u, 

-2U 2 + ‘2UV 1 , 


(O/LO) (61) 


where the only NO correction is in the power spectrum 
itself, 

5A 2 + ^{k) = ^ qi U, (NO). (62) 

Here again OLO has the same accuracy as NO since 
gi(lna;i) = 0. In particular the consistency relation 
r ss —8nt remains unchanged regardless of the evalua¬ 
tion point In Xf. 

Finally, the OLO approximation provides a simple in¬ 
verse relationship between the observables r, n s and a 
and the potential parameters 


U 

Vi 

V 2 


n s - 1 
2 


3r 

16’ 


a 1 3 r 2 

-I— (n s — l)r H- 

2 4 V s ; 64 


(OLO). (63) 


Even if r is not accurately measured, strong upper limits 
where r <C 8|1 — n s |/3 allow a reconstruction of Vi and 

V 2 . 

Equivalently, those observables allow a local recon¬ 
struction of the potential to cubic order around the field 
value where x((f> o) = x\ 


V 

Vo 


1 + Vu(0-M + ^(0-M 2 + ^(0-M 3 - (64) 


The potential amplitude Vg, or the energy scale of infla¬ 
tion, can likewise be determined from measuring A+ x 
or equivalently A^ and r as is well known. 


IV. OPTIMIZED NEXT- AND SECOND-ORDER 
APPROXIMATION 

The OLO approximation of the previous section is ac¬ 
curate at first order in the deviations from scale invari¬ 
ance 0(1/N) and next-to-leading order (NO) in a hier¬ 
archy of slow-roll parameters separated by AN, the tem¬ 
poral scale of features during inflation. For cases where 
AN ~ 1 and the hierarchy is weakly convergent or if 
higher accuracy is desired for future observations, we can 
generalize this approach. In §IV A we consider next-to- 
NO in the hierarchy and an 0(1/N 2 ) in the de Sitter 
deviations. In §IVB, we develop a new optimized NO 
(ONO) approximation that retains the simplicity of the 
NO approximation. This also allows us to contrast OLO 
and ONO with the standard second order (SO) approach 
which conflates AN and N in §IV C. 


In terms of the hierarchy of slow-roll parameters, the 
next-to-NO approximation involves keeping terms to or¬ 
der p = 2 in the Taylor expansion of G in Eq. (17) 


2 

I 0 ~ G(lnxf) + ^2 g P G^(lnx/). (65) 

p= i 


To also include all 0(1/N 2 ) effects, we evaluate the 
second-order terms I\ and J 2 from Eq. (14). In these 
corrections, G' and f'/f can be taken as constants eval¬ 
uated at In Xf as their evolution introduces terms of 
0(1/N 2 AN) which we omit. Combining these pieces, 
we obtain 


lnA^(fc) ~ G(lna;/) + q-^G' (\nxf) + q 2 G" (\a.Xf) 


+ -[G'( lnx f )} 2 -4 


y( ln Af) 


1 2 


( 66 ) 


The same expression is valid for tensors with the appro¬ 
priate replacements as before. 

Since the variation in In k of the last two I\ and J 2 
terms in Eq. (66) is itself 0{\/N 2 AN), we do not in¬ 
clude their impact on n s or a. The only 0(1/N 2 ) terms 
that appear for those quantities come from the nonlinear 
relationship between G64 and slow-roll parameters. This 
justifies our inclusion of such terms in a in the OLO ap¬ 
proximation since there are no further corrections from 
intrinsically second-order GSR effects. 

As shown in the Appendix, in terms of the potential 
slow-roll parameters we obtain 


A 2 


V 


12t x 2 U L 


1 + ( 3gi - - ) U - 2 9l Vi - 2 q 2 V 2 


„ 2 1 03 37t 2 \ 

392 - 3"1 - T + — I" 


+ ( -4 q 2 + -qi + 15 - 2?r 


UV i 


a : 


3 2 13 27r2 , 

1 - m + 2 Vi + 2 qi v 2 + 2 q 2 v 3 

+ (^qi - y j U 2 - ^8<?1 - 0 UV 1 + ^Vi 2 , 
—2V 2 - 2<7iV 3 - 2g 2 V 4 - 6U 2 + 8WVi, 


(67) 
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for the scalar expression and 


A 


2 

+,X 


r 


n t 


at 


V 
67 r 2 


1 + I 9i - g )U 


I 4 8 7T 2 \ 9 

+ |-92 + 3«i-3 + ¥ )« 


+ ( 2<?2 — 2gi + — ) UV 1 


8 U -16 qi U{U - Vi), 

-W + (^r - ^ U 2 - 2( gi - lJWVx, 
-2W(W-Vi), ( 68 ) 


for the tensor expression. These are the master equations 
for all next-to-NO in the hierarchy and second order in 
1/N approximations. They of course also include lower- 
order approximations such as LO, OLO and SO with the 
appropriate zeroing of terms. 


of V 2 back to cfr(lnxi ) from </>(In £ 2 )- We can therefore 
define instead the reconstructed parameters as the V p of 
Eq. (59) which can equally well be predicted from any 
given model for comparison. Moreover in the quadratic 
terms the difference between employing V p and V p is at 
most 0(1/N 2 AN) which we neglect. Thus the master 
equation can be written entirely in terms of V p . Iterat¬ 
ing the inversion we obtain 


U S 

r 3 i/ 

Vi p 

n s — 1 3r 
2 + 16 


+ (— - — 
v 16 96 

v 2 - 

a 1 . 

= -2 + 4 ( ”*- 


9i 2 
64 r ’ 


12 

(n s - l)r + 
3 r 2 

1,r+ 64' 



(ONO), (71) 


for the potential reconstruction. 


B. Optimized Evaluation 


We can again optimize the evaluation epoch In x/ to 
make the optimized NO approximation as accurate as 
a general next-to-NO approximation. We therefore take 
cfe (In £ 2 ) = 0 in Eq. (25) and pick the solution that min¬ 
imizes the next correction |g 3 |, 


<72 (In £ 2 ) = 0 , 


9i (In ^ 2 ) 

hi X‘2 


a/ 37r 2 — 4 

6 


0.84, 


lnxi — 3 i(lnx 2 ) ~ 0 . 22 , 


(ONO). (69) 


We call the evaluation of Eqs. (67) and ( 68 ) with these 
values the optimized next-order (ONO) approximation. 
Note that this has the effect of zeroing the highest V p 
term in each observable. 

For a strong hierarchy AN 1, despite the shift in 
the actual evaluation point to lna: 2 , observables remain 
effectively evaluated at lnaq since 

I 0 = G(lnx 2 ) + gi(lnx 2 )G , (lnx 2 ) + ... 

= G(lnx 2 ) + (lnxi — In z^G'(In x 2 ) + ... 

= G(lnxi) — -(lnxi — lnx 2 ) 2 G ,, (lnxi) + ... 

= G(lnxi) + 32 (In X\)G" (Inxr) + ... (70) 

where ... represents terms and higher. Evaluation 
at lnx 2 simply increases the accuracy by including the 
first correction to this approximation. Since lnxi is a 
maximum of |g 2 | this correction can become important 
for AN ~ 1. Furthermore lnx 2 is a minimum of |< 7 3 1 and 
so even the G^ error is minimized in ONO (see Table I). 

In terms of the reconstruction of G, Vi, V 2 from n s , a , r 
the system is no longer closed due to the appearance of 
31 V 3 in a. As in the case of / 0 above, this term repre¬ 
sents the leading-order effect of a shift in the evaluation 


C. Standard Second Order 


We can also contrast the OLO and ONO approxima¬ 
tions with the standard second-order approach. The 
standard second-order approach is not simply a choice 
of In Xf = 0 in the master equations. By conflating AN 
with N, it truncates at a fixed order in both. In other 
words it assumes V p = 0(1/N P ) and so in the 0(1/N 2 ) 
expansion, 


v 3) v 4 ,... = o, (SO), 

3 ! » 1.06, 32 » 0.21 (72) 


in Eqs. (67) and ( 68 ). Thus, the standard second-order 
(SO) approximation implicitly assumes that the poten¬ 
tial is cubic in form (see the Appendix for the relationship 
to the second-order Hubble flow approximation). Since 
truncation at a fixed V p is not the same as truncation at 
a fixed order in the Taylor expansion, this means that 
different observables are effectively evaluated at an in¬ 
consistent lnx/ or field position. We shall see that this 
makes SO an even worse approximation than OLO for 
cases where AN <C \N\ and the potential is not cubic. 

In particular, dropping 31 V 3 terms in a but keeping 
3 iV 2 terms in n s — 1 leads to the reconstruction 


W «^3i ( 3V, 

8 8 v ’ 64 ’ 


Vi 


V 2 


n s - 1 3r 31 (n s - l ) 2 
2 + 16 + 2 “ 12 

( iw ( 3< 1 1 a_ 7 
16 96/ ; V 12 ^ 768 

a 1/ ^ 3 r2 

^2 + 4 (n ^ 1)r+ ^4’ (S ° } ’ 


(73) 


which is not simply a redefinition of 31 in Eq. (71). 

In particular a now appears in the formula for Vi. For 
a purely cubic potential, these differences simply reflect 
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a different expansion point <f>o for the reconstruction in 
Eq. (64). However for a potential that contains quartic 
corrections, the SO reconstruction can lead to inconsis¬ 
tencies between the parameters {U 1 Vi,X > 2 ) that do not 
simply reflect evaluation at a different field value for the 
potential. We now turn to an illustrative example of such 
a model. 


V. MONODROMY CASE STUDY 


Axion monodromy inflation [17] provides a well- 
motivated theoretical model where the various approx¬ 
imation schemes discussed in the previous sections can 
be tested. The monodromy potential is the sum of a 
smooth and oscillatory part 

V(4>) = V(<I>)+6V(<I>), (74) 

for which we take the simplest case 

v{4>) = x</>, 

SV(4>) = A 4 cos + d^j . (75) 

Depending on the frequency of this oscillation, the cur¬ 
vature power spectrum may gain a large local running of 
the tilt [18, 20] or high frequency features for which the 
slow-roll parameters have an inverted hierarchy and lose 
their ability to characterize observables directly [21]. 

In §V A, we study the low frequency case and test the 
accuracy of the LO, OLO, ONO and SO approximations. 
We develop optimized GSR approximations for the high 
frequency case in §VB, extending results of Ref. [21] for 
second-order corrections in the large amplitude limit as 
well as establishing the connection of the running ampli¬ 
tude, frequency, and phase of the oscillation to parame¬ 
ters of the potential. 


A. Low Frequency Oscillations 


In the low frequency regime, the slow-roll parameters 
are hierarchical and yet the running of tilt can be large 
and nearly constant across the well-measured few efolds 
of the cosmic microwave background (CMB) and large- 
scale structure. In this case, the OLO and ONO slow- 
roll expressions should provide a good approximation for 
observables whereas we have argued that the standard 
SO expressions do not. 

For tests of these approximations, we choose the pa¬ 
rameters of the potential such that the smooth piece pro¬ 
vides a power spectrum amplitude and tilt of the right 
order at some fiducial scale ko, 

A \{k 0 ) = A s = 2.5 x 10~ 9 , 

n s = 0.97, (76) 



-4 -2 0 2 4 


ln(k/k 0 ) 

FIG. 1. Low frequency u) = 1/3 oscillations in the curva¬ 
ture power spectrum A \ under the leading-order (LO) and 
optimized leading-order (OLO) approximations. While the 
LO approximation introduces a large error in the power spec¬ 
trum compared to numerics, the OLO results show that they 
can be corrected by simply shifting the relationship between 
lnfc and (j>. Here the amplitude of the oscillation is set to 

Oi max — 0.01. 


which determines 



(77) 


As shown by the OLO approximation, fluctuations freeze 
out at 


k 0 r)((j > 0 ) = xi « 2.89, (78) 


which provides a conversion between ko and <f >o which we 
use as a definition throughout. The parameter / deter¬ 
mines the frequency of the oscillation and thus the sepa¬ 
ration between terms in the slow-roll hierarchy. Since 


„ J_ _ \dt/dN\ 
Vp~i ffo / 


(79) 


we fix / by a choice of u> ~ AY -1 . Given Ecp (54), 
the series expansion in p converges if uj < 2. Thus in 
the low frequency regime w<C 1, the optimized slow-roll 
approach should be a good approximation. 

We set the phase of the oscillation so that a is at an 
extrema a max at ko, 


0 = ^ - Mod(^ 0 //,27r), 


(80) 
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-4 -2 0 2 4 


ln(k/k 0 ) 

FIG. 2. Low frequency to = 1/3 oscillations in n s , a and r 
under the LO and OLO approximations as in Fig. 1. The 
OLO approximation again provides an excellent approxima¬ 
tion that corrects the shift in the LO results even in this case 
where to, the hierarchical separation between a and n s — 1 , is 
not much less than unity. 

for low frequency oscillations. Finally we set the ampli¬ 
tude of the oscillations A 4 by relating it to the running 
of the tilt at the ko extrema 

A 4 = ~A<^/ 3 a max . (81) 

In practice we choose a max = —0.01 in our examples. 
Since A 4 increases as w -)■ 0 at fixed a max and n s — 1, 
higher-order terms in 5V eventually become important. 
Therefore, we use a max simply as a proxy for choosing A 4 
rather than associating it with the actual value of a at k 0 . 
Likewise, A s and n s also receive higher-order corrections 
and are used here simply to set parameters. 

By solving the exact Mukhanov-Sasaki equation (1) 
with Bunch-Davies initial conditions (5), we can compare 
the power spectra observables to the various slow-roll ap¬ 
proximations and integral approaches. 

We start by considering an example with to = 1/3. 
In Fig. 1, we compare the standard LO and OLO ap¬ 
proximation of Eq. (60) against the numerical calculation 
for the curvature power spectrum A)/. The OLO ap¬ 
proximation reproduces the numerical calculation at the 
10~ 3 level, while the LO approximation misestimates the 
power spectrum by a considerable amount. On the other 
hand, we know analytically that the two differ simply in 
that the power spectrum is shifted by Alnfc ss — In x\. 
Since the relationship between k and <fi depends on the 
physics of reheating, this type of error just introduces an 



FIG. 3. Low frequency to — 1/3 oscillation trajectories in 
the n s -a plane under the LO and OLO approximations as 
in Fig. 1. Both LO and OLO trajectories are good approx¬ 
imations to the numerical results emphasizing that the LO 
results differ simply in the evaluation point along the tra¬ 
jectory, parametrically evaluated here for the same range in 
observed wavenumber k. 


efold shift in those inferences rather than on the shape 
of the potential. 

These considerations apply to n s , a, and r as illus¬ 
trated in Fig. 2. The OLO approximation is highly accu¬ 
rate with differences from the numerics at the 10 -3 level 
or less, whereas the LO results are all shifted by the same 
amount as the power spectrum. 

This joint shift of all observables means that for the LO 
approximation all four observables are consistent with 
arising from the same potential as would be inferred from 
the better OLO approximation. For example, we show 
trajectories in the n s -a plane in Fig. 3. The difference 
between the two approximations now appears simply as 
a different starting and ending point on the same trajec¬ 
tory. 

The importance of consistency becomes clear when 
comparing potential reconstruction between the LO and 
OLO approximations in Eq. (63). In Fig. 4, we show 
reconstruction using the exact, numerical results for n s , 
a, and r as a function of In k as if they were precisely 
measured from data. Again the OLO approximation is 
highly accurate for reconstruction, whereas the LO ap¬ 
proximation shows a shift in tp due to the change in the 
evaluation epoch 

Consistency in the evaluation epoch between observ¬ 
ables leads to consistency in the potential reconstruction. 
In Fig. 5, we show trajectories in the V 1 -V 2 plane. The 
difference between the LO and OLO approximations is 
again the starting and ending point on the trajectory. 
Both approximations would reconstruct potentials that 
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-0.4 -0.2 0 0.2 0.4 


0-00 

FIG. 4. Low frequency w = 1/3 potential reconstruction in 
U, Vi, V 2 under the LO and OLO approximations as in Fig. 1. 
Here the numerical results for n s ,a,r are used in each case. 
The OLO approximation is accurate at the 10~ 3 level whereas 
the LO approximation produces a systematic but consistent 
shift in field position. 

were consistent with the true potential at the same level 
of accuracy. Though the values of Vi and V 2 that they 
recover do differ, they are consistent with a shifted eval¬ 
uation on the potential. 

This should be contrasted with the standard second- 
order (SO) approximation. In Fig. 6, we show the SO and 
ONO approximations for A^ from Eqs. (67) and (68). 
Both approximations do an excellent job of correcting 
residual errors in the LO and OLO approximations in 
the power spectrum itself. For n s ,a,r, Fig. 7 shows that 
ONO further improves the error in their predictions from 
OLO to well below 1C)- 3 . Residual errors are largest for 
a due in large part to the neglect of the 0(1 / N 2 AN) 
term 8 q\UV 2 in Eq. (67) and would decrease for smaller 
choices of to. On the other hand, SO only corrects LO for 
n s — 1 and r leaving the large fractional error in a that 
was present at LO. 

Ironically this means that the SO approximation per¬ 
forms worse than LO in an important sense. In Fig. 8, 
we highlight the problem for the trajectories in the n s -a 
plane. Unlike the LO approximation, the SO approxima¬ 
tion predicts pairs of these observables that are strongly 
inconsistent with the true trajectory. 

Likewise for potential reconstruction, the ONO ap¬ 
proximation improves the reconstruction to better than 
the 10~ 4 level. In fact by defining fully self-consistent 
corrected potential parameters V n from Eq. (59), recon¬ 
struction is actually more accurate than the power spec¬ 



-0.01 0 0.01 


V, 

FIG. 5. Low frequency to = 1/3 oscillation trajectories in 
the V 1 -V 2 plane under the LO and OLO approximations as 
in Fig. 1. Both LO and OLO trajectories are good approx¬ 
imations to the numerical results emphasizing that the LO 
results err simply in the evaluation point along the trajectory, 
parametrically evaluated here for the same range in observed 
wavenumber k. 


trum prediction since it absorbs 0(1/N 2 AN) terms into 
the definition. On the other hand, while the SO approx¬ 
imation corrects the shift of the LO approximation in IA 
and Vi, it does not for V 2 as shown in Fig. 9. The con¬ 
sequence is that in the V 1 -V 2 plane shown in Fig. 10 the 
SO approximation reconstructs points that are not con¬ 
sistent with being anywhere on the trajectory of the true 
potential. This can lead to incorrect falsification of the 
true model from observations. 

Perhaps surprisingly, the OLO and ONO approxima¬ 
tions still perform well for uo = 1 where all terms in the 
infinite V n series are of the same order. This is related 
to the fact that the Taylor series still converges, albeit 
slowly, for even larger values up to w < 2. In Fig. 11, we 
show their predictions for n s — 1, a and r. For the pre¬ 
cision of current measurements, the OLO approximation 
actually still suffices. For higher precision measurements, 
the ONO approximation now makes relatively important 
corrections to the scalar observables and performs well 
at the 10“ 3 level. 

Finally, the accuracy of the OLO and ONO approxima¬ 
tions are directly related to the accuracy with which the 
local tilt and running describe the global power spectrum 
across the observed efolds through Eq. (38). In Fig. 12, 
we show that in the to = 1/3 case where the OLO approx¬ 
imation suffices this approximation is very accurate over 
a large range in efolds. For the to = 1 case, the approxi¬ 
mation only holds for the central ~ 4 efolds and likewise 
the ONO corrections to OLO become important. Thus, 
if the observations are a good fit to the local model, then 
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FIG. 6. Low frequency u> = 1/3 oscillations in the curvature 
power spectrum under the standard second-order (SO) 
and optimized next-order (ONO) approximations. In both 
cases, the second-order corrections are sufficient to correct 
errors in Fig. 1 to below the 10 -3 level (bottom panel). 


FIG. 7. Low frequency to = 1/3 oscillations for the error in n s , 
a and r relative to numerics under the SO and ONO approx¬ 
imations. ONO further improves OLO but SO only corrects 
LO in n 3 and r compared with Fig. 2. a remains shifted in 
SO and breaks the internal consistency of the observables. 


OLO is a good approximation, and ONO provides the 
means to improve its accuracy if necessary in the future. 
If the local model is a bad fit to the data, then AN <1/2 
(or u! > 2 ), and no slow-roll hierarchy of parameters can 
describe the spectrum accurately. 


B. High Frequency Oscillations 


For ui > 2, the slow-roll hierarchy does not converge 
with a finite number of terms. Nonetheless, the GSR 
expansion remains valid so long as the deviations from 
scale invariance, controlled by the amplitude not the fre¬ 
quency of the oscillation, remains small. We must, how¬ 
ever, evaluate the GSR integrals in Eqs. (7) and (14) di¬ 
rectly. Given that in the high frequency limit V p increases 
with p 1 we can no longer naively use their relationships 
with G M given in the Appendix. Reference [14] showed 
that, even in the high frequency limit, 

G' = 3U- 2Vi + € H 0(G') + 0(61). (82) 


We can therefore establish the leading-order connection 
to the potential in the usual way. Namely for the oscil¬ 
latory part, 


6 G = 



9 § v w 

R(i) ’ 


(83) 


is linear in the highest derivative term (or 62 ) and 
hence linear in A 4 in that limit. This extended valid¬ 
ity will be important in the consideration of quadratic 
terms below. It can be demonstrated by integrating 5G' 
assuming that the field rolls according to the smooth po¬ 
tential. Even large amplitude high frequency oscillations 
give highly suppressed effects on the field position which 
we can neglect. 

Under this approximation for the roll, we can convert 
field position to time relative to some fiducial evaluation 
point 77 * as 

<p ~ </>* + V^Hv/v*) ~ </>* + — In (x/kr]*) (84) 

0 * 


to obtain the leading-order contribution from the oscilla¬ 
tions <51nA^(/c) « 6 Io(k), where 


6I 0 



dx 


W'(x) 


sin (a;*In 2 ; + ip) 


+ —— cos (w*lna; + ip) , 

Z0* 


(85) 


with w* 


(/</*) 1 and 
k ) = y 


w* ln(fc77*) + 9. 


( 86 ) 


Moreover, this leading-order result remains valid even As in the optimized slow-roll calculations, the epoch 77 * 
for large amplitude, high frequency oscillations since G' and the associated field value </* can be optimized to 
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FIG. 8. Low frequency w = 1/3 oscillation trajectories in the 
n s -a plane under the SO and ONO. ONO improves OLO but 
SO is worse than the LO (cf. Fig. 3 and also Fig. 14 for the 
second-order Hubble flow result). 


make the approximations as accurate as possible for a 
given k. It does not have to coincide with the fixed ij o 
and <f> o, associated with horizon crossing for fco- 

We can evaluate both terms in the integral using 
Eq. (19), 


dx 

x 


W'(x)e 


xu )* In x 


= - 2 " 


cosh 


7rw*\ 3T(2 + zw*) 




(1 - zw*)(3 - zw*) 


/ 9t coth(7 ru;«/2) i(7r/2 -p) 

V 2(9+ w2) ' 


(87) 



FIG. 9. Low frequency u> = 1/3 potential reconstruction error 
in U, Vi, V 2 under the SO and ONO approximations. ONO re¬ 
constructs the potential to extremely high accuracy (~ 10~ 4 ), 
but SO corrects only U and Vi compared with Fig. 4 leaving 
V 2 with large fractional errors. For ONO, we plot <5Vi, SV 2 - 


Putting these results together, we obtain 


6 I 0 (k) = -A 


cos(-i/> — p) 


3/ 
2 </>* 


sin(^) — j3) 


(90) 


where 


A = 


12 A 4 

A/pi + (3/</>*) 2= 



The phase factor /?(w*) can be expressed exactly in terms 
of the r function formula. Note that it takes on simple 
limiting forms 

lim fj = —-cu* In an, 

uj * —>-0 2 

/3oo = lim = w* [1 — ln(a;*/2)] — y ■ (88 ) 

Here lnaq « 1.06 is the same factor that enters into 
the OLO approximation for reasons that will be clear 
below. The large w* limit corresponds to the resonant 
or stationary phase solution of x = w*/2 for the phase 
factor 2x — w* lnx [21]. To clarify the connection between 
these limits, it is useful to have a simple approximation 
to the exact /3(w*) and we find that 


/3(w*) ~ /3oo + 


7tt 1 

4 yj\ + 1.59w°- 92 + 


(89) 


This should be compared to the result from Ref. [21], 


SA^ = A s 



n s — 1 

A cos [wln(fc/fc 0 ) + v 3 ] > 


(92) 


where <p is considered a free phase parameter. Since this 
form is primarily used for the purpose of fitting a func¬ 
tional form for the deviations to data, in this comparison 
differences in the form are more important than those in 
the association of A, w, and with parameters in the 
potential. On the other hand, our treatment has the ad¬ 
vantage of preserving all of those relations. 

For example, Ref. [21] ignores corrections of order 
//</>* ~ 2 en/to* as they are not relevant for high fre¬ 
quency oscillations, but these may be absorbed into a 
redefinition of amplitude and phase. More importantly 
they take 




to better than 0.03 accuracy everywhere. 


5 In A\ 


(93) 
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FIG. 10. Low frequency ui = 1/3 oscillation trajectories in 
the V 1 -V 2 plane under the SO approximation as in Fig. 1. SO 
is worse than LO (cf. Fig. 5) and can lead to an incorrect 
falsification of the true model. 



which is valid in the |ff.| <C 1 limit but causes differences 
that cannot be reabsorbed into fits at large A. Note 
that Eq. (92) does not guarantee a positive definite power 
spectrum. 

In terms of the phase, Eq. (92) defines only the k de¬ 
pendence rather than the absolute relationship between 
<p and 6 given by our treatment. The dependence of the 
phase on k in Eq. ( 86 ) is consistent with Eq. (92) if we 
choose the evaluation epoch as 77 , = 1 /fco, 

if) m — u;* ln(fc/fe 0 ) + 9. (94) 

The advantage of the more general Eq. (90) is that we 
are not required to do so. Using a fixed evaluation epoch 
omits the evolution of quantities like Eh between that 
epoch and the true freeze-out epoch which depends on k. 
In the more refined treatment of Ref. [22], these effects 
are absorbed into a running or drift of the frequency. 
Note that this type of running does not invalidate our 
use of the slow-roll approximation in Eq. (84) which only 
requires constant Eh during freeze-out of a given k rather 
than across many efolds of k and hence N space. 

We can instead parallel and extend our low frequency 
treatment by allowing the freeze-out epoch for the oscil¬ 
latory piece 6 V to depend on k in potentially a different 
way than for V. We start by reexamining the low fre¬ 
quency case where the two freeze-out epochs should be 
the same. Indeed, if we take the evaluation point as 
fc? 7 * = Xi, the lncci phase shift in /3 in Eq. ( 88 ) cancels 
with the opposite shift in if) in Eq. ( 86 ) so that 



-4-20 2 4 


ln(k/k 0 ) 


FIG. 11. Intermediate frequency ui = 1 oscillation in n s , a, 
and r under the OLO and ONO approximations. OLO is 
still sufficiently accurate for current measurements, but ONO 
provides notable corrections and preserves 10 -3 accuracy in 
observables. 


Equation (90) in the w<Cl limit is then exactly the OLO 
approximation using Eq. (83): 


SIq 


^G\ x — Xl , 
2A 4 


(OLO) 





(96) 


Note that here </>* is a function of k due to the opti¬ 
mized evaluation, and the explicit dependence can either 
be directly computed given the background solution or 
characterized by the running of the frequency from some 
central value at ko [ 22 ]. 

For high frequency oscillations w > 1, we increase ac¬ 
curacy by shifting the evaluation point 77 * to the station¬ 
ary phase point 


lim (^ - /?) = ^ - w* + 8 + j, kr]* = L ^~. (97) 

co* —>00 j 4 z 

Optimized evaluation again means that is given by 

the solution to 


W*) = y 


1 

2/V 


(98) 


To cover both low and high frequency cases, we can com¬ 
bine the two freeze-out criteria as 


lim (ip — j3) 


f +9 + 


7r 
2 ’ 


kr 7 * = X\. 


(95) 


fc 77 * = max (aq, -y) . 


(99) 
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FIG. 12. Local, n 3 — a parameterization of the power spec¬ 
trum from Eq. (38) compared to numerics for the low and 
intermediate frequency cases. Accuracy of the local approxi¬ 
mation is directly related to the hierarchical structure of the 
slow-roll parameters and hence that of the OLO truncation 
of the hierarchy. For u) = 1/3 both suffice, whereas for u> = 1 
both require correction from higher terms in their respective 
Taylor series. 


FIG. 13. High frequency to = 20 oscillation in the curvature 
power spectrum under the analytic GSR integral approx¬ 
imation as compared with numerical calculations for a model 
with A ~ 0.2 amplitude oscillations (top). Differences are 
less than 10~ 2 with a significant correction coming from the 
quadratic terms ( Q) which mainly change the phase and am¬ 
plitude of the oscillation as well as shift the normalization 
parameter InA,, (bottom). 


For the second-order corrections in the u> -C 1 limit, 
the ONO approximation applies and captures all terms. 
In the opposite limit ofu> 1, quadratic terms in A can 
be more relevant since a large amplitude oscillation does 
not produce large changes to the well-constrained average 
tilt and running. In this large amplitude and frequency 
limit, the dominant second-order correction is due to I\ 
[14] where a similar resonance occurs in the stationary 
phase approximation. Thus, the oscillatory part of the 
power spectrum is well approximated by 

SI 0 + 2l 1 Sh+SI'l, (GSR) (100) 

where the smooth and oscillatory parts of I\ are given by 

h = ^(l-n s ), ( 101 ) 

and 

1 f 00 dx 

5Ii = — / —SG' (In x)X(x). (102) 

V 2 Jo x 

Repeating the same steps as for /q, we have integrals 


that can be evaluated using 

r°° 

/ — X'(x)e iu * lnx 

Jo x 


= —2 _l “*« sinh 


m 


3r(2 -}- ica*) 

(1 — lw*)( 3 — ioj *) 


/9w*7rtanh(7ru;*/2) 

2(9 + W 2) 6 ’ 


(103) 


to obtain 


A 

dli = tanh 

y/2 



sin(i/> — (3) + 


2 </>* 


cos(t/> — /3) 


(104) 

Note that the Slf correction contains squared or 2w* fre¬ 
quency contributions that cannot fully be absorbed into 
the form of Eq. (92). However to quadratic order in A, 
the combination of the Iq and I\ differences in form be¬ 
tween Eqs. (100) and (92) can be absorbed into a redef¬ 
inition of A Sl A and the oscillation phase for w 1. 
To see this note that, aside from the phase "0^/3, these 
quadratic differences are given by 


Q 


sjl 

2 


+ 2IiSIi + dli 


(105) 


and \d2dli and dlo are related by a phase shift of 
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7r/2 such that their squares sum to a constant since 
tanh(7rw*/2) —► 1. 

As an example, we choose u> = 20 and A 4 = 1.78 x 
10~ 13 which gives A(ko) ~ 0.2. Other parameters are 
set as in the low frequency case. We use the ONO ap¬ 
proximation to determine A s , n s , and a for the smooth 
potential and add its contribution using Eq. (38) to com¬ 
pare the total 


In A n = In 


ONO 


S In A\ 


(106) 


to numerical results. Note that, by implementing the 
ONO approximation, where V is evaluated at kri(<j> 2 ) = 
X 2 , we use two held positions <j>* and <fr 2 for each k which 
allows us to simultaneously optimize for the oscillatory 
and smooth parts of the potential. 

In Fig. 13, we show that even for this large ampli¬ 
tude oscillation, Eq. (100) is accurate to better than 
10~ 2 without any adjustment to the phase, frequency 
or amplitude of oscillations. Note that, as expected, the 
quadratic terms in Q , while significant in establishing the 
good agreement with numerics, takes the form of a con¬ 
stant plus an oscillatory piece of the same frequency as 
the leading-order contribution. Thus, it can be absorbed 
into a redefinition of parameters in Eq. (92) and the de¬ 
pendence of the oscillation phase ip — /3 on k. For even 
larger values of the amplitude A , this is no longer possi¬ 
ble and so accurate expressions will require calculations 
to cubic and higher order in A 


VI. DISCUSSION 

By utilizing the GSR approximation, we have provided 
a systematic study of the evaluation and interpretation 
of power spectra where a relatively large local running 
arises from features in the potential. This approach as¬ 
sumes only that the average deviation from scale invari¬ 
ance is small and associated with the number of efolds to 
the end of inflation as \/N. The frequency of temporal 
features 1/AJV can be much larger than 1/N. We intro¬ 
duced the GSR slow-roll hierarchy parameters G&) which 
are directly related to observables and elucidate their re¬ 
lation to the standard hierarchies of parameters in the 
potential and Hubble flow. This parameterization works 
for P(X, <j>) models as well with a suitable generalization 
of the conformal time to the sound horizon. 

For models with 1 < AN -C N, a slow-roll hierarchy 
of parameters still exists, but the running of the tilt a is 
only smaller than n s —1 by a factor of AN not N. In these 
models, the leading-order slow-roll calculation gains rela¬ 
tively large corrections since a large running implies that 
the slow-roll parameters are not constant. Instead the 
calculation and interpretation of such models proceeds 
by Taylor expanding the temporal evolution. This series 
converges rapidly as long as AN > few, and this con¬ 
dition also guarantees that observable power spectra are 
well characterized by the local tilt and running across 


the well-measured few efolds of the CMB and large scale 
structure. 

In fact next-to-leading-order corrections in the Taylor 
series of each observable can be consistently reabsorbed 
into an optimized temporal evaluation of leading-order 
terms. This is because the first Taylor correction just 
represents a shift in the epoch of fluctuation freeze-out 
to about an efold before horizon crossing. The advantage 
of this OLO approach is that the analysis of observables 
proceeds exactly in the same way as a standard leading- 
order analysis and only the interpretation in terms of 
their correspondence to the inflationary model differs. 
For the interpretation of current data fits to constant 
tilt and running of the tilt, this simple approach suffices 
in accuracy. In terms of potential reconstruction it re¬ 
covers a local cubic expansion evaluated consistently at 
the optimized field point. 

Contrast this with a common approach in the literature 
that attempts to correct for the evolution with a second- 
order expansion under the assumption that AN ~ N. 
In this case the tilt is corrected for evolution between 
horizon crossing and freeze-out but the running of the 
tilt, which is assumed to be intrinsically second order, 
is not. In this case, the observables are effectively eval¬ 
uated at inconsistent epochs. Potential reconstruction 
from a second-order approach likewise produces incon¬ 
sistencies which could potentially lead to an incorrect 
falsification of the true model from the observations. For 
models with large running it can only be consistently 
applied to a purely, rather than locally, cubic potential. 
Ironically then the second-order analysis provides a more 
complicated but less general approach compared with the 
leading-order analysis as we explicitly demonstrate using 
a low frequency axion monodromy model. 

For even higher accuracy, we can keep both 1/iV 2 terms 
and optimized next-to-leading-order terms in the hierar¬ 
chy. This ONO approach leads to better than 10~ 3 accu¬ 
racy in observables and potential reconstruction even for 
AN ~ 1 where all of the infinite hierarchy of parameters 
have the same magnitude. 

Finally for AN <C 1, the slow-roll parameters possess 
an inverted hierarchy and lose their utility for predicting 
observables. Even in this case the GSR approach pro¬ 
vides accurate predictions as long as the amplitude of 
deviations from time translation invariance is small re¬ 
gardless of the frequency of temporal features. Likewise 
the technique of optimizing the epoch of evaluation can 
be generalized to establish precise relationships between 
the potential parameters and observables. 

We test this GSR approach with axion monodromy 
models in the high frequency limit and provide expres¬ 
sions that are accurate to second order in the amplitude 
and optimized in the temporal evaluation. The optimized 
approach also provides the direct relationship between 
parameters of the potential and phenomenological tem¬ 
plate fits to the running phase, frequency and amplitude 
of oscillations in the power spectrum used in the litera¬ 
ture. For sufficiently large amplitude deviations and/or 
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future high precision measurements, relating power spec¬ 
trum features to potential parameters requires going be¬ 
yond the second order calculations presented here. We 
leave such considerations and their impact on interpret¬ 
ing observables to a future work. 
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Appendix A: Parameter Relations 


To relate the three different parameterizations of the 
slow-roll hierarchy G {eh, 8 p } and {U, V p } used in the 
main text, we establish here the relationship between 
the GSR, Hubble slow-roll, and potential slow-roll ap¬ 
proaches. 

We expand expressions to 0(1/N 2 ) under the assump¬ 
tion that 

{'»•") = °(v), 

{G<-Up,V p }=o(-^ij =T ). (Al) 

with 1 <C A N < |AT|. Unlike the standard slow- 
roll approximation, we therefore keep terms that are 
0(1/N AN P ) but still drop those that are 0(1/N 2 AN P ) 
for p > 1. For example we keep V 2 and IAV \ but omit 

uv 2 . 

The relationship between the GSR and Hubble slow- 
roll variables involves the conversion between conformal 
time and efold 


dN 
din 77 


—aHr). 


(A2) 


This quantity can be parameterized by using Eqs. (4) and 
(40) to obtain 


dH 

dN 

d 2 H 

dAP 


—£hH, 


—+ 2Si)H. 


(A3) 


Integrating the Taylor expansion of H around N to find 
V = Jn dN/e N H , we obtain 


(A4) 


Using this relation, the definition of / from Eq. (3), and 
the Hubble slow-roll hierarchy equation (40), we obtain 
for the GSR scalar parameters 


In / 2 

f 

f 


In + 2eH + + 4e J?di, 

—<5i — 2 ch - 4 e 2 h - 3enSi, 


, or 2 r 32 2 28 2 2 

G « 4e h + 2<5i + -82 + — e# + 1 p e u 8 i - -dj, 

G" ~ — 2 A 2 — — 83 — 8 ejj — lOe/fdi + 2d 2 , 

G ( p ) « (_i)p+i ^2d p + ^d p+1 ) , (p > 2) (A5) 

and for the tensor parameters 


ln fh 

fk 

fh 

G' h 

G'L 

n (p) 



+ 2 + 5 e~ h + 4e/fdi, 


-eh - 3e 2 h — 2eh8i, 

n 22 2 16 
% e H + ~^ e H + ~7^£h8i, 

-4 e 2 h — Aeh8\, 

0 , (p> 2 ). 


(A 6 ) 


For the relationship with the potential parameters, we 
iteratively solve Eq. (48) 

U = 2e h + 0(N~ 2 ), 

Vi =e h -8 1 - 6 ^ + 0(N~ 2 ) (A7) 

along with 

V P = - d ^+0(N~ 2 ) 

= (-l)P^8 p + 6 ^j+0(N- 2 ) (A 8 ) 


for p > 2. Notice that this is the same combination that 
enters into G^ in Eq. (A5). Inverting these relations we 
obtain 

« = 0+O(iV- 2 ), 

*i = !-£(5)"'W+0(W- 2 ), 

00 , S n 

Sp = E (3 ) V «+2 + 0(N~ 2 ). (A9) 


We then plug these back into the exact Eqs. (47) and 
(48) and use 


V' dU 

V dcj) 
V'dVp 

V d(f> 


2W(Vi -W), 


aHjj ~ 1 eh T 3e~h T 2eh8\. 


Vp+i + [(p-l)Vi-jW]V p , (A 10 ) 
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to obtain 


become 


3 H 2 

~v~ 

£h 

U 

51 

5 2 

5 P 


U U 2 UV 1 

1 -1 - 1 -— 

6 12 9 ’ 


1 U 

2 ~~ ¥ 

U 


E 


v n + -U 2 - -UV, + 

9 6 


5 yl 

"18"’ 


n—0 


T-E(i) + 


n—0 


EU) K+2+^-Ev-i+V^ 


(-i)'EU v, 


n—0 


/ n+p? 


(P>2) 


(All) 


where we have performed one more iteration of eh to ob¬ 
tain eh/U to 0{N~ 2 ) since it appears in the denominator 
of A 

As noted above, an advantage of the potential slow-roll 
parameters is that in the Taylor expansion for G' mixed 
p terms in S p go away when rewritten in V p , 


Gs 


-V? 

9 

G 's 

17 9 5 2 9 

a 3G - 2Vi + — U 2 - -UV, - -V 2 , 

5 3 3 


G" - 

a -2V 2 -6U 2 + 8UV,, 


G(p) s 

a — 2 V P , (p > 2 ), 

(A12) 


for the scalars and 


G h ? 



G' h * 

a U + |w 2 - 2WVi, 


u h n 

a -2U 2 + 2UV,, 


r>(p) - 

a 0 , (p> 2 ), 

(A13) 


for the tensors. 

With these relations, the full 0(N~ 2 ) expressions with 
AN -C 1 for the curvature power spectrum observables, 


In A\(k) ~ G(lna;/) + ^ q p G^ (In x/) 
p =1 



[G'(lnx f )] 2 


-4 


'f 

./ 


2 


(In x f ) 


n s- 1 ~ -G'{\n Xf ) Yl q p G^ + 1 \lnx f ), 

P= 1 


a sa G"(lna;j) + ^ q p G^ p+ 2 \lnxf), (A14) 

p=i 


A. 


1/ 


n ~ 12tt 2 G . 


1+ (3 gi --)G-2^ gp V p 

p=l 


, 2 103 3tt 2 

+ [ 3, 2 - — 


w 2 


+ ( -492 + ^1 + 15 - 2tt 2 ) UV, 
. 2 13 2-7T 2 1 2 

4,2 - 3 ?1 - y + ^> v , 


1 - 3W + 2V, + 2 E 9pVp+i 

P=1 

17 


a : 


+ - y) u ' ( v 8,?1 ~ 3 ) UVl + 3 Vl ’ 

00 

- 2 V 2 - 2 ^ QpVp +2 - 6 U 2 + mV,. (A15) 

p=i 


The analogous relations for the tensor observables give 


A 2 

^+,x 


4 8 7T 2 

~ 92 + 3 91 _ 3 + ~6 


U 


+ 2^2 — 2gi + 


17 


WVi 


r & 8 U - 16q,U(U - V,), 
nt ~ —G + E - G 2 - 2( gi - 1)GV 1; 


at ~ -2U{U- Vi). 


(A16) 


These expressions can be used to derive any second-order 
approximation specified consistently by the order in the 
hierarchy of q p (lnxf) corrections or in the standard ap¬ 
proach by keeping only G(l/N AN) terms with In a;/ = 0. 

For example, in terms of the Hubble flow parameters 
of Eq. (46) , keeping the conversions in the standard ap¬ 
proach where we assume e n = (D(l/N) gives 


n s — 1 — — 2 ei — e 2 — 2 e 2 + ( 2 q, —— ) £162 


£ 263 , 


+ Ui - 


a = — 2 e,e 2 — E 2 E 3 , 


n t = — 2 e, - 2 e\ 
a t = - 2 e,€ 2 , 


2 < 7 i - ) eie 2 , 


(HSO) 


(A17) 


which reproduces the standard result used in the litera¬ 
ture (e.g. [1, 6 , 7]) withqi = gi(0) = 7/3 — In 2 — 7 ^. Note 
that because the truncation is inconsistent in the models 
we consider with AN <C \N\, these relations differ nu¬ 
merically from the SO approximation expressed in terms 
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of the same truncation in the potential parameters in the 
text unless the potential is cubic. Nonetheless n s — 1 and 
a are still inconsistently evaluated, and the Hubble flow 
second order (HSO) also performs worse than the simpler 
OLO approximation (see Fig. 14). 


FIG. 14. Low frequency to = 1/3 oscillation trajectories in 
the n s -a plane under the OLO and HSO approximations as 
in Fig. 8. HSO like SO provides an inconsistent trajectory 
and is worse than the simpler OLO approximation. 
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